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ABSTRACT 

We present the analysis of an approximately 3 year long Rossi X-ray Timing Ex- 
plorer (RXTE) monitoring campaign of the canonical soft state black hole candidates 
LMC X-l and LMC X-3. In agreement with previous observations, we find that the 
spectra of both sources can be well-described by the sum of a multi-temperature disk 
blackbody and a power law. In contrast to LMC X-l, which does not exhibit any pe- 
riodic spectral changes, we find that LMC X-3 exhibits strong spectral variability on 
time scales of days to weeks. The variability pattern observed with the RXTE All Sky 
Monitor reveals that the variabili ty is m ore complicated than the 99 d or 198 d period- 



icity discussed by Cowley et al. ( |199l| ). For typical ASM count rates, the luminosity 
variations of LMC X-3 are due to changes of the phenomenological disk blackbody 
temperature, fcTm, between ~ IkeV to ~ 1.2 keV. During episodes of especially low 
luminosity (ASM count rates ^ 0.6 counts sec -1 ; four such periods are discussed here), 
fcXln strongly decreases until the disk component is undetectable, and the power law 
significantly hardens to a photon index of F ~ 1.8. These changes are consistent with 
state changes of LMC X-3 from the soft state to the canonical hard state of galactic 
black holes. We argue that the long term variability of LMC X-3 might be due to a 



wind-driven limit cycle, such as that discussed by Shields et al. (1986) 



Key words: accretion - black hole physics - Stars: binaries - X-rays: Stars - Stars: 
LMC X-l, LMC X-3 



1 INTRODUCTION 

Long term variability on time scales of months to years is 
seen in many galactic black hole candidates. In analogy to 
the 35 d cycle of Her X-l, the long term variability of some 
objects has been identified with the precession of a warped 
accretion disk. Possible driving mechanisms for a warp in- 
clude radiation pressure due to the luminous central X-ray 



source ( Pringlc 1996[ Maloncy et al. 1996 ; Maloney & Begel- 
man 1997] ; |Maloney et al. 1998|), torques exerted by an ac- 



cretion disk wind (Schandl & Meyer 1994; schandl 1996), or 



tidal forces (Larwood 1998, and references therein). 



In Isome sources, the long term changes in the X-rav 
luminosity have been associated with state changes of the 
accretion disk. At low luminosities, these sources are usually 



observed in the hard state, in which the X-ray spectrum is 
dominated by a hard power-law component with a photon 
index of T ~ 1.7 and an exponential rollover at ~ 150 keV. 
The hard state spectrum is usually described in terms of 
thermal Comptonization (Sunyaev & Triimper 1979; Dove 
et al. 1997; and references therein). At higher luminosities, 
black hole candidates exhibit a soft spectrum that can be 
characterized by a (multi-temperature) blackbody with a 
peak temperature of kT ~ 1 keV. In addition, a power law 
with a photon index of F ~ 2.5 or softer is present. Changes 
between the spectral states are typical for black holes as is 



evidenced by the 1996 soft state of Cygnus X-l (Cui et al 
1997; Cui et al. 1998), and the frequent transitions seen 



in GX 339-4 (Makishima et al. 198rj; Bouchet et al. 199; 



Wilms et al. 1999a; Corbel et al. 2000). See Tanaka & Lewin 
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(1995), Nowak ( 1995), and references therein for further de- 



tails. 

Prior to the observations discussed here, LMC X-l and 
LMC X-3 were the only black hole candidates that had al- 
ways been seen in the soft state. Thus, both were ideal can- 
didates for a systematic study of the properties of the soft 
state. Although both sources were otherwise thought to be 
spectrally very similar, LMC X-l had not shown signs of any 
periodic long term variability, whil e LMC X-3 was known to 
be vari able on a ~100 d time scale ( |Cowley et al. 1991 ; ^ 



ley et a . 1994 ). We therefore initiated a two resp. three year 
long campaign to observe LMC X-l and LMC X-3 with the 
Rossi X-ray Timing Explorer (RXTE) . During the campaign 
we performed X-ray observations of roughly lOksec length 
at approximately three week intervals. The spacing and ex- 
posure time of the observations were chosen such that they 
would enable us to track spectral changes of the sources 
over any long term variability pattern. The campaign was 
started in 1996 December with two long (~ 170ksec) obser- 
vations of the sources. The analysis of these observations is 
presented in a companion paper (Nowak et al., 2000, hence- 
forth paper I) . Preliminary results from the first year of this 
campaign, using earlier versions of the response mat rix and 
backgr ound models, have been p resented elsewhere ( Wilms 



et al. 1999b; Wilms et al. 1999c), this paper is devoted to a 



discussion of the first two years of the campaign. 

The remainder of this paper is structured as follows. 
In Section |^ we present the details of our data analysis 
procedure. Section |§] is devoted to the study of LMC X-3. 
Section ^| contrasts these observations with those from 
LMC X-l. In Section ^ we interpret our results in the con- 
text of current models for the soft state of galactic black 
hole candidates, and in the context for models of the long 
term variability of X-ray binaries. We summarize our results 
in Section ^. Throughout this paper we assume a distance 
to the LMC of 50kpc. 



final analysis was obtained from these individual matrices 
by adding the matrices weighted by the fraction of photons 
coming from each of the intervals. We analyzed data taken 
in the energy band from 2.5 keV to 20keV. 

Since our data are background dominated above 
~5keV, good background modeling is essential for our anal- 
ysis. Background subtraction of the PCA data was per- 
formed using a model taking into account sky pointings 
of the PCA, and modeling the background variability us- 
ing data from the Very Large Event (VLE) counter of the 
PCA or using the so-called "Faint model" , depending on the 
source count rate. To estimate the uncertainty of the back- 
ground model we compared the count rate measured in the 
PCA channels above 30keV, where no source photons are 
detected, with the background model flux at these energies. 
In all cases the agreement between the background model 
and the measured background was within 2%. We accounted 
for this fluctuation by renormalizing the background model 
flux, using the xspec corrfile facility, such that x 2 was min- 
imized in our spectral fits. Note that this introduces a "hid- 
den" fitting parameter to the data analysis. This parameter 
was not taken into account in the number of degrees of free- 
dom since we assume that it converges to the correct nor- 
malization of the background model, and therefore does not 
introduce an uncertainty into the spectral fitting process. 

We used data from the ASM in order to be able to 
place our observations in the context of the long term vari- 
ability of the sources. This instrument, an array of shadow 
cameras scanning the whole sky visible from the spacecraft 
for five to ten tim es per day (Remillard & Levine 1997; 



Levine et al. 1996), provides almost uninterrupted informa- 



tion about the long term behaviour of bright X-ray sources. 
We used the "definitive one-dwell" data available from the 
Goddard Space Flight Center (GSFC) and only used flux 
solutions for which Xrcd < IT. 



2 DATA ANALYSIS 

Onboard RXTE are two pointed instruments — the Pro- 
portional Counter Array (PCA) and the High Energy X- 
ray Timing Experiment (HEXTE) — as well as the All Sky 
Monitor (ASM). We used the standard RXTE data analysis 
software, ftools 4.2, to examine the PCA and HEXTE data. 
Spectral mode ling was done using xspec, version lO.OOab 
(Arnaud 1996). Due to the short duration of the pointed 
observations we use only the PCA data for this analysis. 
We used essentially the same data screening and analysis 
strategy as in paper I, i.e., we used PCA top xenon layer 
data only and ignored data taken within 30 minutes after 
passages through the south Atlantic anomaly and where the 
background count rate, as measured by the "electron ratio" , 
was comparatively large (see paper I for details). Data from 
all PCUs was combined for the final analysis. Contrary to 
paper I, the observations presented here are so short that 
the Poisson error dominates the uncertainty of the spec- 
trum; therefore, no systematic error was applied to the data. 
For some of the observations, only part of the PCA detec- 
tors were turned on. These intervals were extracted sepa- 
rately and then combined. Response matrices were gener- 
ated for each of these intervals. The response matrix for the 



3 LMC X-3: LONG TIME SCALE SPECTRAL 
VARIABILITY 

3.1 Introduction 

The black hole candidate LMC X-3 was discovered during 
UHURU ob servations of the Large Magellanic Cloud (Leong 
et al. 1971). A summary of the early observational history 
has been given by Treves et al. ( 1988 ) . Interest in LMC X-3 
was heightened when it was realized that its X-ray spectrum 
is very si milar t o that of Cyg X-l in the high state. White & 
Marshall ( 1984 ) therefore established the source as a poten- 



tial black hole candidate. Subsequent radial veloci ty mea- 
surements, resulting in a mass function of 2.3 Mq( Cowley 
et al. 1983), as well as optical photometry (van der Klis et 
al. 1983; Kuiper et al. 1988; and references therein) of the 
B3V compa nion lead to a minimum mass of ~4MQ( Mazeb 
et al. 1986), wi t h a most probab le mass of 9-10 Mq (Cow- 



[Ley et a~ 198E ; Paczyhski 1983 ) for the compact object 
LMC X-3 thus almost certainly contains a black hole. 

The luminosity of LMC X-3 can be as much as ~30% of 
the Eddington luminosity of a 9 Mq black hole. The source 
has been known to be strongly variable since Ariel V and 
the High Energy Astrophysics Observatory 1 (HE AO 1) dis- 
covered large (> 5) intensity changes on time scales of weeks 
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to months ( |Griffiths fc Seward 1977| ; |johnston et al. 1979| ). 
A first systematic study of its Io ns term periodicity was per- 
formed by Cowley et al. (1991), who, based on Ginga and 
HEAO 1 data, found a long term period of 98.9 d or 197.8 d. 
In addition, optical photometry indicated a lag of ~ 20 d 
between the optical and the X-rays (the optical leading the 
X-rays; Cowley et al. 1991). Later Hubble Space Telescope 
(HST) investigations re vealed the presence of the periodicity 
also in the ultraviolet (Cowley et al. 1994), albeit the long 
term trends were not as clear as in the earlier data. These 
HST observatio ns an d the shape of the X-ray lightcurve led 
Cowley et al. (1994) to prefer the ~ 99 d period and to 
postulate that LMC X-3 might have gone through a low- 
luminosity state sim ilar to the "extend ed low" states that 
are seen in Her X-l (Parmar et al. 1985). 

Detailed spectral studies of the long term variability 
were first performed by Ebisawa et al. ( 1993] ), who used 
data from Ginga LAC pointings, taken over a period of three 
years, that covered a large range of source luminosities. Ebi- 
sawa et al. showed that the X-ray spectra could be well- 
described with a multi-temperature disk blackbody model 
plus a power-law spectrum. Alternatively, optically thick 
(r c ?t 20) Comptonization could be used (see also Treves 
et al. 1988), but is not preferred due to the unphysical pa- 
rameters resulting from the spectral fitti ng. In accordance 
with Cowley et al. ( |l99l| ), Ebisawa et al. ( |l993[ ) found that 
the soft spectral component hardens with increasing source 
luminosity (i.e., its temperature inceases), while the power- 
law tail was found to be independent of the luminosity. One 
possible interpretation of this result is a varying mass accre- 
tion rate M. 

A major shortcoming of the earlier spectral observa- 
tions is that it was impossible to put the observations in 
the context of the overall long term variability of the source 
due to the lack of a real all-sky monitoring instrument. Fur- 

t hprmore thp rplativp inflWihilit.y of thp parly X-ray sat.pl- 
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lites m ade it impossible to obtain an even sampling of X-rav 
spectra over the long time scale changes observed. A satellite 
such as RXTE, which combines a broad spectral coverage 
with very flexible scheduling and an all-sky monitoring in- 
strument, is therefore ideal to systematically study the long 
term behaviour of sources such as LMC X-3. In the remain- 
der of this secti on we discuss the RXTE-ASM soft X-ray 
li ghtc urve (§3.2), describe our analysis of the po inte d data 
(§3.J), and present the results of our campaign (§3.4). 



3.2 Long Term Variability of LMC X-3 

In Fig. |lj: we display the long term light curve of the object 
as observed with the ASM, binned to a resolution of 3d. 
Note the strong variability on time scales of ~ 100 and ~ 
250 d, which is well sampled by our monitoring (see Fig. |^). 
This variability had already been noted during the first year 
of ASM operations (Levine et al. 1996). 

The variability pattern itself is very complicated. The 
analysis of the AS M ligh t curve using the g eneralized peri- 
odogram of Lomb (1976) and Scargle (1982) reveals signifi- 
cant periodicities on long time scales of ~ 100, ~ 160, ~ 190, 
and ~ 240 days (Fig. | l|b). Epoch folding analysis (Leahy 
Schwarzenberg-Czerny 1989| ; Davies 1990) inde 



et al. 1983 




Figure 1. a) Upper left: Dynamical PSD of LMC X-3 for the 
ASM data from 1996 through mid 1999. The white line dis- 
plays the behaviour of the average ASM count-rate (a period of 
350 d corresponds to 4 ASM cps). b) Upper Right: Lomb-Scargle- 
Pcriodogram for the ASM data through 1999, the dashed line 
indicates a false alarm probability for periods at the 0.01% level 
(for the definition of the false alarm probability, see Scargle 1982). 
c) Lower left: ASM lightcurve of LMC X-3, binned to a temporal 
resolution of 3 d and variation of the spectral parameters for our 
RXTE pointed observations. See text for explanations. 



by these methods on the 1.7 d orbital time scale of the sys 



tern. This is consistent with earlier results (Weisskopf et al 



1983; Cowley et al. 1983; van Paradijs et al. 1987), but con 



pendently verifies this result. No significant variability (us- 
ing the full temporal resolution ASM lightcurve) is detected 



trary to other high mass X-ray binaries, such as Cyg X-l or 
Vela X-l, where a clear orbital modulation of the X-ray flux 
is detected in the ASM. For these latter sources, this vari- 
ation is most pr obably due to p hotoelectric absorption in 
the stellar wind ( Wen et al. 1999 ). The lack of orbital vari- 
ability in LMC X-3 can therefore be interpreted as a lack of 
a strong stellar wind from the companion star. This is not 
very surprising, considering that the companion is a B3 V 
star dWarren fc Penfold 19751 ; |Cowley et al. 1994]) ■ 

Compared to sources with clear perioditicies such as 
Her X-l, the peaks in the Lomb-Scargle periodogram of 
Fig. Qb are quite broad. This can be a sign that the periods 
detected by this method are only quasi-periodic, i.e., they 
change with time or are not present during the whole time 
span covered. We therefore computed a dynamical Lomb- 
Scargle periodogram by taking slices of the ASM light curve 
(binned to a resolution of 3d) with a length of 730 d each, 
and shifting these slices with a step size of 9 d over the whole 
ASM data rangeQ For each of these slices the Lomb-Scargle 
periodogram was computed for the period range from 6 d to 



* In order to extend the dynamical periodogram to the full length 
of time covered by the ASM data, we added gaussianly distributed 
random data spanning one year before and after the ASM light 
curve. Each fake data set was chosen to have the same mean 
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365 d. To ensure that the periodograms are comparable, we 
normalized them such that 



Table 1. Observing log of the monitoring observations of LMC 
X-3. 



a 
/< 



PSD(/) df 



(1) 



where / m i n and / max are the minimum and maximum fre- 
quency for which the Lomb-Scargle periodogram is com- 
puted, and where a 2 is the variance and /i the mean ASM 
count rate of the 730 d lightcurve. The resulting periodogram 
was then gray scale coded and is displayed in Fig. [l]a, with 
each periodogram being displayed at the mid-time of the 
light curve for which it was computed. 

Although the individual Lomb-Scargle periodograms 
are obviously not statistically independent, our approach is 
useful in revealing the long term trends of the light curve. 
Fig. ^a shows clearly the origin of the individual peaks in 
the total periodogram of the source. Albeit with varying sig- 
nificance, the dynamical periodogram shows that a ~ 100 d 
periodicity is present in all data segments analyzed through 
the middle of 1998, confirming the periodicity discussed by 
Cowley et al. (1994). The ~ 100 d periodicity is mainly at- 
tributable to the times of low (< 0.5 ASM cps) source lumi- 
nosity, that are, e.g., seen in 1996 July, 1997 January, and 
1997 August (Fig. |c). After mid 1998, the amplitude of the 
100 d periodicity is seen to decrease; however, some of this 
reduction in amplitude can be attributed to the increasing 
inclusion of fake data in a given two year periodogram. A 
possible shift in the distance between the light curve minima 
during mid 1998 through 1999 may lead to the additional 
time scale of ~ 160 d present in Fig. Qb (although see the 
discussion below). The ASM lightcurve for 1996 looks simi- 
lar to the folde d HE AP 1 and Ginga lightcurve presented by 
Cowley et al. (1991). This might be an indication that the 
"quasi-sinusoidial" variability of 1996 is not just a random 
event, but occurs frequently. 

The peaks at lower frequency in the total Lomb-Scargle 
periodogram are due to long term periodicities that are not 
fixed in period. The dynamical Lomb-Scargle periodogram 
shows a dominant periodicity on a time scale of ~ 190 d in 
1996, which increases to ~ 270 d for the 730 d light curve 
centered on 1997 December. These shifts in period can be 
attributed to the maxima of the ASM light curve, which 
behave quasi-sinusoidally in 1996 and then show an increase 
in period in 1997 and the first half of 1998, before (possibly) 

returning to a more periodic behaviour^]. The strength of the 
individual peaks in the total periodogram (Fig. ^b) can thus 
be seen as partly depending on the duration of the intervals 
in which each period dominated. It is also possible that the 
~ 160 d peak seen in the periodogram is the continuation 
of the ~ 190 d period seen at the start of 1996, while the 
~ 240 d peak is due to the formation of a new periodicity 
beginning in late 1996 or early 1997. It is also interesting to 



and variance of the first and last year, respectively, of the actual 
measured data set. 

T Note that the determination of periods <; 250 d in the dynam- 
ical Lomb-Scargle periodogram is difficult since it is based on 
segments of 730 d length only. Since the long period shows up in 
the total periodogram of Fig. [l]b, however, we believe the dynam- 
ical periodogram to represent at least the general trends in the 
period shifting behaviour. 



Obs. 


Date 


Exposure 


Count Rate 




ymd 


sec 


counts s — 1 


01 


1996.12.05 


139800 


390.1 ±0.1 


02 


1996.12.29 


7800 


94.3 ±0.1 


03 


1997.01.16 


7600 


83.6 ±0.1 


04 


1997.02.07 


9300 


123.7 ±0.1 


05 


1997.03.09 


8000 


352.2 ±0.2 


06 


1997.03.21 


7200 


374.6 ±0.2 


07 


1997.04.17 


8500 


318.5 ±0.2 


08 


1997.05.05 


5200 


319.3 ±0.3 


09 


1997.05.27 


8400 


289.4 ±0.2 


10 


1997.06.18 


8100 


246.0 ± 0.2 


11 


1997.07.04 


6200 


234.8 ± 0.2 


12 


1997.08.03 


5900 


148.7 ±0.2 


13 


1997.08.19 


8100 


35.9 ±0.1 


14 


1997.09.09 


9400 


102.6 ±0.1 


15 


1997.09.12 


7900 


113.3 ±0.2 


16 


1997.09.19 


8500 


113.4 ±0.1 


17 


1997.09.23 


10300 


153.7 ±0.1 


18 


1997.10.11 


7300 


127.7 ±0.2 


19 


1997.11.02 


8400 


152.3 ±0.2 


20 


1997.11.23 


8900 


111.5 ±0.1 


21 


1997.12.12 


3300 


166.8 ±0.3 


22 


1998.01.06 


8200 


268.0 ±0.2 


23 


1998.01.24 


6100 


254.0 ± 0.2 


24 


1998.02.20 


9100 


213.4 ±0.2 


25 


1998.03.12 


6600 


182.4 ±0.2 


26 


1998.04.06 


9200 


190.0 ±0.2 


27 


1998.05.07 


7600 


169.0 ±0.2 


28 


1998.05.29 


9400 


9.1 ±0.1 


29 


1998.06.29 


3400 


42.5 ±0.2 


30 


1998.07.20 


9800 


137.3 ±0.1 


31 


1998.08.12 


7700 


133.6 ±0.1 


32 


1998.09.02 


4200 


264.7 ±0.3 


33 


1998.09.30 


7600 


33.2 ±0.1 



Exposure times shown are rounded to the closest 100 sec. The 
count rate is background subtracted 



note that the ~ 100 d and ~ 190 d periodicities that exist at 
the beginning of the ASM light curve have approximately 
the same period separation as the ~ 160 and ~ 240 day 
periods that exist at the end of the ASM light curve. 

Although the ASM does in principle provide X-ray 
colour information, LMC X-3 is too weak to enable us to 
detect any statistically significant changes in the spectral 
shape during the three years of ASM coverage. Usable X- 
ray colours can only be computed after rebinning the data 
on ^ 4d intervals. These data indicate that the source 
might harden during lower ASM count rates. With the ASM, 
however, it is not possible to characterize these changes in 
greater detail. We therefore turn to our pointed monitoring 
observations to further describe these spectral trends. 
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Figure 2. Typical fit of a multi-tcmpcrature disk spectrum and 
a power law to data from LMC X-3. The bottom panel shows 
the residuals in units of a. The residuals are consistent with the 
remaining calibration uncertainty of the PCA. 



3.3 Spectral Modeling of LMC X-3 PCA Data 

Using the screening criteria of section ^, we extracted the 
data from our 1997 and 1998 observations. A log of the ob- 
servations is given in Table Our conservative screening 
criteria reduce the effective usable exposure times from the 
scheduled lOksec to the times shown in the Table. 

We modeled the RXTE spe ctra using the standard 



multi-t emperature disk model (Mitsuda et al. 1984 



Mak- 



ishima it al. 1986) plus an additional power law component. 
This is the traditional spectral model for describing the X- 
ray continuum of the soft state. We further discuss below 
the interpretation and possible shortcomings of this model 
(§P, see also paper I). Our fits generally gave acceptable 
results. Several observations, however, showed evidence for 
deviations of the residuals in the iron Ka band. Although 
these deviations are mainly in the 2a range, adding a narrow 
(<r = 0.1 keV) iron line to the model improved our fits. We 
therefore added a narrow Gaussian line feature at 6.4 keV to 
all of our fits. Although the presence of a line is consistent 
with our analysis from paper I, we note that Fe line stud- 
ies with the PCA are complicated by the pre sence of sys- 



tematic features in the PCA response matrix (Wilms et al 
1999a|^|and by the fact that the power-law and the disk 
model have comparable flux in the Fe-line region. The line 
parameters that we found here, therefore, should be seen as 
consistent with an upper limit of ~60 eV for the equivalent 
width. We took the absorption in the intervening interstellar 
medium into account by fixing the equivalent hydrogen col- 
umn to the value found from radio and ASCA observations 
N-a = 3.2 x 10 20 cm~ 2 (Lister-Staveley, 1999, priv. comm.; 
see also paper I), usingthe cross sections of Baluciriska- 



vations had the large Ah values that are detectable with the 
PCA energy range and spectral resolution (Ah ^ 10 22 cm -2 ; 
Stelzer et al. 1999). 

We list the results of our spectral modeling in Tab. g, 
and a typical spectrum is shown in Fig. [| Generally, the 
spectral model described the data satisfactorily with re- 
duced x 2 values of Xred ^ 2. We note that especially for some 
brighter observations, the formal \ 2 values are influenced by 
systematic features in the response matrix. This can be ex- 
amined, e.g., by comparing the ratio between the data and 
the best fit model, with the ratio between an observation of 
the Crab pulsar and nebula and the best fit power law to 
that observation (Wilms et al. 1999a). Visual inspection re- 



Church & McCammon (1992). Using a fixed Ah is justified 
since preliminary modeling showed that none of the obser- 



veals that the deviations between the data and the model for 
all monitoring observations were comparable to or smaller 
than those seen in the Crab. Thus, although Xicd > 1 in 
these cases, adding additional spectral components to the 
model would only serve to model the response matrix fea- 
tures and not the observed spectrum. We are therefore confi- 
dent that, given the spectral resolution of the PCA and the 
faint flux level of LMC X-3, we have found a satisfactory 
spectral description for all pointed observations. 

We emphasize that the interpretation of the spectral 
components in terms of physical parameters, however, is 
hampered by the properties of the multi-temperature disk 
blackbody model. This discussion is necessitated by the re- 
sult that the typical power law indices found in our data 
analysis are softer than those found by Ginga (F ~ 2.2; 
Ebisawa et al. 1993). A more detailed discussion can be 
found in paper I. A problem associated with using the multi- 
temperature disk blackbody plus a power law for modeling 
the soft state is that especially for large V (i.e., soft power- 
law tail) , the power law can have a non-negligible flux in the 
lowest pulse height analyzer (PHA) channels so that the fit 
parameters of the power law are almost completely deter- 
mined by these low channels. This effect is present in most 
of our fits for which r <; 4. We tried to limit the influence 
of the power law in the higher PHA channels with several 
other spectral models, e.g., by adding a heavily absorbed 
power-law to the disk blackbody, by using a modified multi- 
temperature disk blackbody which includes a high energy 
power law tail only above a given energy, or by applying 
thermal Comptonization models. None of these models re- 
sulted in acceptable x 2 values. In addition, using the multi- 
temperature disk blackbody and leaving the equivalent col- 
umn Ah a free fit parameter always resulted in Ah — > cm 2 . 

The only model, apart from the one that we adopted 
here, which allowed us to describe the spectrum of LMC X-3 
with similar or even better x 2 values was a multi- 
temperature disk blackbody to which a low energy black- 
body spectrum was added. Such a model has been success- 
fully applied, e.g. , also to the 1997 RXT E observations of 
GRO J1655-40 ( [Rothschild et al. 1999b - Although such a 
model works for those observations of LMC X-3 where we 
find F ^ 4, the model fails to reproduce those observations 
where a harder power law component is clearly present in 
the data. We therefore decided to use the multi-temperature 
disk blackbody spectrum plus a power law in our analysis. 

The significant contribution of the power law in the low- 
energy PHA channels has some effect on the determination 
of the characteristic disk temperature, fcTi n . As we noted in 
paper I, the blackbody temperature obtained from ASCA 
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Figure 3. Temporal variability of the spectral parameters of LMC X-3. 



data was systematically lower than the kTi n derived from 
the PCA data. We speculated that this was due to the sig- 
nificance of the power law in the low energy PHA channels 
and the differing low energy cutoffs between the two detec- 
tors. We further speculated that the ASCA fit temperature 
was more characteristic of the temperature for the "seed 
photons" that are Compton upscattered to become the high 
energy spectrum. As we further discuss below for our obser- 
vations of LMC X-l (jfej), the presence of the power law 



may greatly affect the determination of the blackbody nor- 
malization, y4disk. 



3.4 Transitions between the Soft and Hard State 

We now discuss the variations of the spectral parameters 
in terms of the long term variability of LMC X-3. In Fig. ^ 
we display the spectral variation of LMC X-3 in the con- 
text of its ASM long term light curve. The plot reveals a 
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Table 2. Results of Spectral Fitting to the LMC X-3 Data. 



Obs. 




kT ltl 
keV 




4-disk 




r 


A PL 

1 n— 1 
1U 


Alpine 
10 — 5 


EW 

cV 


x 2 


/dof 


01 


1 


23 ± 0.01 


33 


9 ± 0.7 


2 


4 ±0.2 


0.7+0-, 

—0.2 


8.5+iV 

— 8.5 


11 


59 


5/47 




n 


09 j- n on 

.JZ ZE U.UU 


7 "1 


q+0.8 
y -l.l 


2 


in u.z 


n c; -I- n 9 


8 4 i 2 2 


six 


79 
t z 


7 /A9 
1 1 *±z 


03 





91 ± 0.01 


29 


7+1-6 
' -1.7 


3 


± 0.1 


1.4 ± 0.3 


6 7+ 2 - 

u. 1 _ 2 7 


65 


60 


8/42 


04 


1 


01 ± 0.01 


29 


n+0.5 
—0.8 


4 


3 ± 0.3 


—1.0 


6.8 ± 2.4 


51 


19 


3/41 


05 


1 


91 +0.00 

Z1 -0.01 


29 


n+0.6 
a -0.7 


3 


1 ± 0.1 


3 4+ - 5 

J - 4 -0.4 


11 6+ 4 - 1 


22 


62 


6/41 


06 


1 


22+0-0? 
— 0.01 


31 


—0.6 


3 


2 ± 0.1 


3.3 ± 0.5 


20.5+ 4 -! 

—4.5 


36 


59 


5/41 


07 


1 


90+0.01 

Z,3 -0.00 


27 


5 ± 0.4 


1 


1 ± 0.2 


5 0+1.5 
— l.i 


7 K+4.1 
'■°-3.4 


16 


7.3 


3/41 


08 


1 


22 ± 0.01 


29 


5 ± 0.5 


3 


7 ±0.2 


0+1.2 

—0.9 


12.0 +5 A ° 

—4.7 


25 


16 


4/41 


09 


1 


20+°°° 
-0.01 


29 


± 0.4 


4 


1 ± 0.2 


4 8+1-4 
-1.1 


13 7 + 2'5 

J 1 —3.7 


32 


19 


9/41 


10 


1 


I 7 +0.01 

II -0.00 


28 


+0.4 
U -0.5 


4 


9+0.3 

z -0.2 


5 5+ 1 - 8 


5 8+ 37 

o.o_ 3 2 


17 


58 


3/41 


11 


1 


16 ± 0.01 


26 


9 ± 0.5 


4 


l+oo 

—0.2 


5 3+ 1 -? 

J,J -1.3 
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31 


58 


1/41 
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6+? 1 
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8 2+ 2 ? 
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°-0.3 


4 6+ 2 - 4 
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71 
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29 


4 ± 0.9 


4 
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3 0+,V? 

—0.7 


6.7+ 2 - 4 


56 


24 


1/41 


16 


1 
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2 1 


7+1.4 
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4 


6 ± 0.3 


6 9+ 29 
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35 


70 
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"1 7 
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UD It U.U1 


97 
z / 
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1 
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U it U.z 
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1 1 
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2 /A 1 
0/ 41 


18 


1 


01 ± 0.01 


29 


6+?- 2 , 

— 1.0 


1 


9+0.3 
— 0.2 


4.1+J-? 

— 1.1 


6.4+ 2 '6 

" — 3.0 


16 


26 


5/41 


19 


1 


06 ± 0.01 


28 


"—0.7 


4 


4 +0.3 
— 0.2 


5.5+J-? 

u — 1.3 


7.5+ 2 -?, 
1 — 3.0 


12 


12 


6/41 


20 


o 


99 ± 0.01 


27 


8+0-5 
-0.8 


4 


°-0.2 


45 +i.8 
1.1 


8.0+ 2 ' 2 


70 


17 


2/41 


21 


1 


08 + 01 


28 


0+2.3 

°-0.8 


4 


O+0.5 
°-0.3 


4 n+3.3 
^■ y -1.5 


11 9+ 47 


59 
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1 
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97 

Z 1 


f\ + a 

O It U.^t 
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7 4- n 9 
It u.z 


7 4+I.9 
'• 4 -1.4 


11 7 4_ q c 
±1.1 It o.D 


1 


1 7 
10 




23 


1 


18 ± 0.01 


27 


3 ± 0.5 


4 


2 ± 0.2 


' —1.6 


11.2+3-? 
—4.5 


30 


38 


4/41 


24 


1 


-0.00 


28 


1 ± 0.4 


4 


4+0-o 
-0.2 


7 8+?? 
' •°— 1.6 


10 6+2? 
— 2.9 


38 


12 


3/41 


25 


1 


10 ± 0.01 


27 


0+0.3 
— 0.6 


4 


g+0.4 
-0.3 


u — 2.0 


10.4 ± 3.4 


15 


35 


2/41 


26 


1 


11 ± 0.01 


27 


k+0.5 
-0.6 


4 


3 ± 0.2 


6.5+H 
u - u — 1.4 


9 3+? n 
^ — 2.9 


38 


55 


8/41 


27 


1 


05 ± 0.01 


25 


o+!i 
u — 1.2 


3 


5 ± 0.1 


4 4+0-1 
^ -0.6 


5.8 ± 3.1 


28 


101 


5/41 


28 










1 


8 ±0.0 


0.1 ± 0.0 


1 2+ 1 - 1 


60 


26 


5/43 


29 





c; fi +0.05 


80 


1+80.5 
i -30.8 


2 


8 ±0.1 


i-i±S:2 


9 n+ 2 - 8 

z ' u -2.0 


33 


22 


7/41 


30 


1 


04 ±0.01 


27 


0±0.6 


4 


tr+0.3 
°-0.2 


6.3t 2 ;0 


4 7+2-6 
— 2.5 


30 


38 


1/41 


31 


1 


16 ±0.00 


24 


g+0.6 
°-0.5 


5 


c+0.8 

°-0.5 


97 7+I8. 1' 

^'■'-10. 5 


7 n+3.9 
' ' u -3.2 


24 


87 


0/41 


32 


1 


17 ±0.01 


27 


9+O.8 
-0.9 


3 


8 ±0.1 


7 7+ 1 - 4 
'■'-1.2 


11 4+ 5 - 2 


30 


43 


3/41 


33 





gg+0.02 
°"-0.03 


57 


n+10.7 
U -7.9 


2 


7 ±0.1 


1.0 ±0.2 


4 n+2.3 
*- y -2.6 


71 


30 


3/41 



Tin, ^4disk : peak multi-temperature disk temperature and normalization. V: photon index of the power law Apl: Power law 
normalization (photons keV -1 cm -2 s — 1 at 1 keV). j4Li ne : Line normalization (photons cm -2 s — 1 in the line), the Gaussian line was 
fixed at 6.4 keV with a width <r of 0.1 keV. EW: line equivalent width. Uncertainties are at the 90% confidence level for one interesting 
parameter (Ax 2 = 2.71), the interstellar equivalent column was fixed at TVjj = 3.2 X 10 20 cm -2 . 



clear correlation between the total source luminosity, as ex- 
pressed by the ASM count rate, and the parameters of the 
soft spectral component. During times of higher source lumi- 
nosity, the temperature of the multi-temperature disk black- 
body, fcTin, increases. At the same time, the normalization 
of the disk blackbody stays remarkably constant (see also 
Fig.|d). This effect is similar to that found in the soft-state 



of other galactic black holes ( Tanaka fc Lewin 1995 ). Dur- 



ing episodes of comparably low ASM count rate, however, 
the spectrum shows clear deviations from the typical pure 
soft state behaviour. During our campaign, there were four 
of these events: in 1997 January, 1997 August, 1998 June, 
and 1998 October. During these times fcTi n decreased to 
fcTm 0.5 keV, while the photon index hardened (Fig. |5|b). 



For the most extreme case, Obs28, no evidence for the soft 
component is seen, indicating that fcTj n 0.5 keV, and the 
photon index is T — 1.8. These values are remarkably similar 
to the typical hard state behaviour seen, e.g., in Cyg X-l. 
We interpret the decreases in the ASM countrate as evidence 
for transitions from the normal soft state into the hard state 
of LMC X-3. 

In Fig. ^ we display the unfolded photon spectrum of 
Obs28, together with the spectra of the three subsequent 
monitoring observations, showing the transition from the 
pure hard state spectrum back to the soft state. Obs28 is the 
only observation of our campaign in which there is no de- 
tectable evidence for the soft component. The occurance of a 
pure hard state is thus either a very rare event for LMC X-3, 
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Figure 4. Spectral evolution of LMC X-3 from the hard state of 
Observation 28 to the normal soft state seen in Observations 30 
and 3f. Shown are the unfolded photon spectra, multiplied by 
photon energy, the lines denote the best fit model. 



Figure 6. The disk blackbody temperature, kT[ n as a function of 
the 1 day averaged ASM count rate. Spectra observed when the 
ASM count rate was less than about 0.6 counts sec -1 (indicated 
by the dashed line) exhibit signs for a transition to the hard state. 



or, if it is associated with each of the dips observed by the 
ASM, it is of relatively short duration. The data presented 
here represent the first observational detection of a pure 
hard state in LMC X-3. We speculate, however, that the 
data taken on HE AP 1 day 461 and mentioned by White 
& Marshall ( 1984 ) might have been observed during one of 
the soft to hard state transitions. During this observation, 
however, the source was too faint for HEAO 1 to obtain 
even a reliable measurement of the hardness ratio. Further- 



more, during 2 000 April another hard state was seen (Boyd 
& S maje 2000|), with th e timing behavior being hard state 
like ( [Homan et al. 200C| ). 

The hard state observation Obs28 is the observation 
with the lowest PCA count rate of our campaign, it was 
also performed close to one of the lowest ASM fluxes seen 
so far. In Fig. |^ we display kTi n as a function of the ASM 
countrate averaged over 1 d worth of ASM dwellings centered 
on the observation. In the figure we have identified spectra 
for which fc71 n < 0.8 keV. All of these spectra also show a 
power law component which is relatively strong compared 
to the soft spectral component. These data indicate that for 
ASM count rates ^ 0.6 cps a transition to the hard state 
behaviour is probable. Therefore, it is possible to use the 
ASM count rate as a trigger for future systematic studies of 
the soft to hard transition in LMC X-3. 

Our monitoring observations are too short to be able to 
characterize the variability of LMC X-3. We computed the 
average rms variability of the source from the background 
subtracted 16 s lightcurves. Although there is an apparent 
trend of the rms variability to increase with lower source 
flux, this trend is fully consistent with systematic variations 
due to PCA background model fluctuations at the 2% level. 



This level is consistent with the uncertainty of the back- 
ground model found in the time series analysis presented in 
paper I, and also with the adjustments to the background 
flux found during the spectral modeling process. It is there- 
fore impossible to check our in terpretation of Obs28 being a 
hard state or an intermediate ( Mendez fc van der Klis 1997 ) 
state from its temporal characteristics. Much longer hard 
state observations (> lOOksec) are needed for this determi- 
nation. 

Outside of the hard states, LMC X-3 exhibits the stan- 
dard behaviour of a classical soft state source. Fig. ^ displays 
several correlations between the spectral fit parameters. For 
power law indices F ^ 4 there is a clear correlation between 
the power law flux and T: the power law component hard- 
ens with decreasing power law fl ux (Fig. Sa). Such behaviour 



has been noted for GX 339- 4 flWilms et al. 1999a|) and the 
Seyfert 1 galaxy NGC 5548 flC hiang et al. 200CQ , albeit in 
their hard spectral states. Figs. ||b and c verify the results of 
Ebisawa et al. ( 1993] ) that during the soft state there is no 
strong correlation between the soft and the hard component. 
For r 4, r does not correlate with the accretion disk tem- 
perature, fcTm. There is a slight trend for the power law flux 
and F to increase with kT[ n ; however, the statistics of the 
data are not good enough to claim any firm correlation be- 
tween these two parameters. In conjunction with the above 
spectral results for LMC X-3, in paper I we noted that the 
Fourier frequency-dependent coherence function indicated 
that the soft and hard X-ray components were uncorrelated 
with each other for LMC X-l (the variability of LMC X-3 
was too weak to perform such an analysis). 

Finally, we consider the behaviour of the soft spec- 
tral component. The normalization of the multi-temperature 
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Figure 5. Correlation between the spectral parameters for LMC X-3. The dashed line in subfigure d denotes the proportionality 
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disk blackbody, Adisk, is nearly constant throughout the 
monitoring observations (Fig. ||d), while fcTj n strongly 
varies. A similar behaviour has been found in those soft state 
black hole candidate s where the luminosity is dominated by 
the soft component (Tanaka & Lewin 1995). If j4disk is inter- 
preted literally, then A^isk oc (ri n /d) 2 cos 9, where r- m is the 
innermost radius of the accretion disk, d is the distance of 
the source, and 9 is the disk inclination angle. In Fig. |5]d we 
display the values of ri n cos 1//2 9 corresponding to our mod- 
eled values of ^diak on the upper a;-axis. The parameters 
for ri n found here are typi cal for galactic black hole candi- 
dates ( Yaqoob et al. 1993j ). Interpreted in this framework, 
the constancy of A^isk for LMC X-3 could imply that the 
geometrical configuration of the accretion disk does not ap- 
preciably change during the soft state episodes of the object, 
and that all changes in the disk luminosity are accounted for 
via temperature changes. We also note that for disk black- 
body temperatures kT[ n 0.9 keV, the data approximately 
agree with the proportionality of Adisk ex T-~ 8 ' 3 which is 



expected if Tin oc r ; ~ 3 ^ 4 and the inner disk radius moves 
outward, with relatively little bolometric luminosity change, 
during a soft-to-hard state transition (Fig. dashed line). 

The literal interpr etation of A^isk and kTj n , however, is 
not without problems (Merloni et al. 2000). The innermost 



region of thin accretion disks is more complicated than as- 
sumed in the theory underlying the multi-temperature disk 
blackbody model. Even in simple a accretion disk models 
( [ghakura fc Sunyaev 1973 ) , the inner region of the accretion 
disk is dominated by electron scattering. Thus the emerging 
spectrum is slightly Comptonized. As a result, the observed 
kTi n does not correspond to the temperature of the inner 
edge of the disk, but rather equals the colour temperature 



of the disk (Ebisawa et al. 1991; Shimura & Takahara 1995) 



Correcting for this effect is stron gly model dependen t, see , 
e.g., Shimura & Takahara ( |l995|) and Merloni et al. ( |200C| ) 
for discussions. Furthermore, the presence of a hard spectral 
component can also influence A«sk, as we mentioned above 
(see also paper I and Sect. ^^). We therefore consider it dan- 
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gerous to use correlations between a (corrected) kT[ n and ri n 
fit parameters to draw firm conclusions on the geometry or 
physical environment of the observed system. The apparent 
independency of the disk normalization from its character- 
istic (colour) temperature kT[ n , however, is striking in the 
case of LMC X-3. Any accretion disk theory attempting to 
model this object will have to explain this independency. 

Finally, we note that no clear correlation is found be- 
tween the equivalent width of the iron line and the other 
spectral parameters. This is consistent with our interpreta- 
tion of the iron line parameters as upper limits. 



4 LMC X-l: LONG TIME-SCALE SPECTRAL 
VARIABILITY 

4.1 Introduction 

Like LMC X-3, LMC X-l was also discovered during the 
UHURU scans of the LMC. Its optical companion is a 0(7- 
9)111 star with a mass function of 0.144 Mq and an orbital 
period of 4.2 d (Cowley et al., 1995, Hutchings et al. 1987, 
Hutchings, Crampton & Cowley, 1983) in a photoionized He 
II nebula (Bianchi & Pakull 1985; Pakull & Angebault 1986) 
The luminosity of LMC X-l is about 2 x 10 3S er| 



Table 3. Observing log of the monitoring observations of LMC 
X-l. 



et al. 1981 ) and was found to be quite constant (pyunyaevl 



et al. 199C) 



The X-ray spectrum of LMC X-l is similar to that of 
LMC X-3 (Ebisawa, Mitsuda & Inoue, 1989; paper I and 
references therein), although the relative flux of the disk 
black body with respect to the power law is smaller. There 
is again evidence for a weak Fe emission line. 

The long term behavior of LMC X-l during 1996 as seen 
by RXTE has been analyzed by Schmidtke, Ponder & Cow- 
ley ( 199!: ). In agreement with previous observations these 
authors did not find evidence for any systematic variability 
on long timescales. We will, therefore, not re-analyse these 
data but concentrate in what follows on our monitoring data 
from 1997 and 1998. 



Obs. 


Date 


Exposure 


Count Rate 




ymd 


sec 


counts s" 1 


01 


1996.12.06 


128900 


151.6 ±0.0 


02 


1996.12.30 


8600 


146.7 ±0.2 


03 


1997.01.18 


8200 


180.4 ± 0.2 


04 


1997.02.08 


9300 


149.4 ±0.1 


05 


1997.03.09 


3800 


157.0 ±0.2 


06 


1997.03.21 


8200 


160.8 ±0.2 


07 


1997.04.16 


10600 


163.4 ±0.1 


08 


1997.05.07 


8600 


163.6 ±0.2 


09 


1997.05.28 


7000 


147.2 ±0.2 


10 


1997.06.18 


3600 


161.9 ±0.2 


11 


1997.07.04 


3400 


148.6 ±0.3 


12 


1997.08.01 


8500 


157.1 ±0.2 


13 


1997.08.20 


7600 


171.1 ±0.2 


14 


1997.09.09 


9600 


138.9 ±0.1 


15 


1997.09.12 


8300 


163.3 ±0.2 


16 


1997.09.19 


9100 


120.8 ±0.1 


17 


1997.10.10 


9100 


154.3 ±0.2 


18 


1997.11.01 


8800 


157.0 ±0.2 


19 


1997.11.23 


6400 


156.1 ±0.2 


20 


1997.12.12 


5100 


160.7 ±0.2 


21 


1998.01.04 


8700 


189.4 ±0.2 


22 


1998.01.25 


7100 


168.0 ±0.2 


2.3 


1998.02.20 


9200 


177.1 ±0.2 


24 


1998.03.12 


3300 


173.9 ±0.3 


25 


1998.04.07 


9400 


178.1 ±0.2 


26 


1998.05.06 


5500 


154.0 ±0.2 


27 


1998.05.28 


8900 


179.1 ±0.2 


28 


1998.06.28 


3700 


178.1 ±0.3 


29 


1998.07.19 


9300 


151.4 ±0.2 


30 


1998.08.13 


8000 


181.2 ±0.2 


31 


1998.09.02 


5200 


143.9 ±0.2 


32 


1998.09.29 


8600 


112.7 ±0.1 



4.2 Long Term Variability of LMC X-l 

The data of the monitoring observations of LMC X-l were 
treated in the same way as those for LMC X-3. A log of 
the observations can be found in Tab. ^. During the cam- 
paign, the count rate of LMC X-l was on average a factor of 
~ 2 smaller than that of LMC X-3. Again, we modeled the 
PCA data using a multi-temperature disk blackbody plus 
a power law and an iron line. We assumed an equivalent 
hydrogen column of Ah = 7.2 x 10 21 cm -2 (Lister-Staveley, 
1999, priv. comm., see also paper I), as is appropriate from 
the radio measurements. The results of the spectral model- 
ing are shown in Tab. ^. 

As for LMC X-3, the LMC X-l spectra are well- 
described by a phenomenological model consisting of a disk 
blackbody (with kT in » 0.8-1 keV) and a soft (r > 3.5) 
power law component. Fig. |^ displays the correlations among 
the fit parameters for LMC X-l. The most striking aspect 
of Fig. |7J is the apparent strong correlation between the phe- 
nomenological disk blackbody temperature, kT[ n , and the 
disk blackbody normalization, Adisk- Furthermore, there is 
a weaker correlation between the disk temperature and the 



Exposure times shown are rounded to the closest 100 sec. The 
count rate is background subtracted 



2-20 keV power law flux. These correlations are such that 
the disk temperature decreases by ~ 20% while its normal- 
ization increases by a factor of ~ 4. This yields an overall 
increase of the blackbody flux by a factor of ~ 1.6. There is 
at the same time, however, a decrease by a factor of ~ 2 in 
the power law flux. In contrast to LMC X-3 the contribu- 
tion of the power law makes up a greater fraction of the total 
2-20 keV flux in LMC X-l; therefore, the decrease [increase] 
in disk temperature [normalization] corresponds to an overall 
decrease in 2-20 keV flux by only » 20%. 

There are two possibilities to explain these spectral cor- 
relations. First, it is possible that, as in LMC X-3, the overall 
decrease in 2-20 keV flux corresponds to the beginnings of a 
transition to a low/hard state. Arguing against this notion, 
however, is that we do not see any corresponding hardening 
of the power law flux. Second, a more likely scenario is that 
we are seeing a complicated systematic interaction between 
the phenomenological disk component and power law com- 
ponent. For the soft photon indices seen in LMC X-l there 
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Table 4. Results of Spectral Fitting to the LMC X-l Data. 



Obs. 


keV 


^disk 


r 


A PL 
10" 1 


^Line 
KT 5 


EW 
eV 


X 2 /dof 


01 


89 + 00 


63 8 + 


3.4 ± 0.0 


1.1 ± 0.0 


10.7 ± 0.0 


63 


80.8/41 


02 


0.96 ± 0.01 


33.0+2'J 


3.4 ± 0.0 


5.9 ± 0.6 


10.8tf° 


53 


84.9/42 


03 


n m + n 1)1 

L/.f iJ _1_ \J.\J1. 


45 2+ 2 - 6 


„ -, +0.0 

• 3 - i -o.i 


4.5 ± 0.5 


19 2+ 3 '° 

ly '^-3.1 


85 


88 2 /42 


01 


0.87 ± 0.01 


63 8 +3 / 2 

—4.1 


3.5 ± 0.1 


7 8+J- 2 

' -1.0 


13.1+ 3 '^ 


67 


50.3/41 


05 


91 + 01 


47 i+ 3 - 5 


3.7 ± 0.1 


8.9 ± 1.8 


8 1+ 4 ' 3 


17 


37.5/41 


06 


0.88 ± 0.01 


61.3+ 2 S 

u — 3.7 


3.6 ± 0.1 


6.7 ± 1.1 


8.6+ 3 'J, 


51 


45.5/41 


07 


0.86 ± 0.01 


69 6+ 2 '5 


3.6 ± 0.1 


7 5+1- 2 
1 1.0 


11 6+2? 


71 


95.8/41 


08 


n as + n ni 

U.OO _1_ U.W-L 


62 5+ 1 ' 7 


3.6 ±0.1 


74 +i.i 
1.2 


11 3+ 2 ' 7 


65 


48.7/41 


09 


0.87 ± 0.01 


55 2+ 3 i 


3.7 ± 0.1 


8.2 ± 1.4 


ii 7+2'2 

±x - 1 —2.9 


76 


52.7/41 


10 


0.85 ± 0.01 


67 5+S i 

u ' ,u — 6.2 


3.7 ± 0.1 


8 9+i 9 

° — 2.0 


13 1+1-2 


78 


32.4/41 


11 


0.83 ± 0.01 


77 6+?; 9 

1 ' '"—11.9 


3.7 ± 0.1 


7 6+ 2 -i 

1 • u -2.0 


6.6 ± 4.2 


16 


37.6/41 


12 


93 + 01 


36 8+ 2 ' 1 


3.5 ± 0.1 


7 7+I.O 
'•'-0.9 


10 2+ 2,8 
lu - -3.0 


56 


91.2/41 


13 


f) ni + f) 01 


53 8 +2 ' 7 


3.7 ±0.1 


8 8+ 1 - 2 
1.3 


q 3+3-1 
y.o_ 3 2 


50 


35.4/41 


14 


0.95 ± 0.01 


29 8+H 
£a -2.4 


3.6 ± 0.1 


8 6+i'S 

°' u — 0.9 


9 4+ 2 ' 9 
a - -2.8 


51 


54.3/41 


15 


0.85 ±0.01 


74.3+?- 9 

^—5.3 


3.8 ± 0.1 


9.0+^ 

— 1.3 


8.6+? 9 

— 2.6 


55 


85.9/41 


16 


89 + 01 


46 5+ 2 ' 3 


3.7 ± 0.1 


g +1.3 
— 1.1 


11. 5+ 28 


73 


68.3/41 


17 


i no 4- o oi 

l.UU _1_ U.W-L 


no i +1.3 
■"- x -1.8 


q c+0.0 


8.0 ± 0.8 


10 5+ 30 


54 


34.5 /4l 


IS 


0.97 ± 0.01 


30.2+i? 

— 2.1 


3.5 ± 0.1 


7.3+°- 9 

— 0.8 


n.o+2 T 

— 3.1 


58 


29.9/41 


19 


n no _i_ n m 

U.yo □! U . U-L 


26 9+ 19 
zo.y_ 2 4 


3.6 ± 0.1 


9.3 ± 1.1 


1 9 c;+ 3 - 3 
LZ -°-3.7 


63 


47.0 /41 


20 


88 4- 01 


53 7+ 3 - 5 
uo - ' -4.5 


3.9 ± 0.1 


10 4+ 2 '° 
lu -1.8 


7 n+3.4 
1 -3.5 


45 


36.9/41 


21 


0.96 ± 0.01 


38 1+ 2 ,'?, 


3.7 ± 0.0 


ii 8+;- 2 
—l.i 


13 1+ 2 ? 

±lJ,± — 3.2 


59 


66.0/41 


22 


0.98 ± 0.01 


29 6+i? 

^"- /,w — 2.2 


3.7 ± 0.1 


9.9 ± 1.1 


10 9+o? 

3.4 


53 


38.6/41 


23 


97 4- 01 


28 0+ 1 ' 8 


3.6 ± 0.0 


10 4+°- 9 

lu -0.8 


14.5 ± 3.0 


66 


90.9/41 


9 1 


n on _i_ n m 
u.yu it U.U1 


c 9 tr + 4.7 
dz,J -8.2 


o. i ± U.l 


y.y ± i.y 


7 7+4.5 
'■'-5.2 


i n 
1U 


on i //ii 
oy . 1 / 41 


25 


0.93 ± 0.01 


37 7+?,'2 

J ' ■ ' -3.0 


3.7 ± 0.0 


11.7 ± 1.1 


10 7+ 3 ,' 2 , 


52 


40.9/41 


26 


0.87 ± 0.01 


58 l+5'I 


3.8 ± 0.1 


9 6+;-^ 


11 8+?'5 
±Xi °— 3.6 


73 


46.4/41 


27 


0.90 ±0.01 


46.8t 3 3 ;6 


q 7+0.0 

■^-'-o.i 


11.8 ± 1.2 


15 9+ 2 ' 9 


79 


44.2/41 


28 


90+ 01 


48.J+" 


3.7 ±0.1 


11.8 ± 1.9 


14 9+ 4 ' 5 


75 


59.5/41 


29 


0.87 ±0.01 


55.8+. 2 ; 9 


3.9 ±0.1 


11 9+ 1 ' 5 


8 1+ 2 ' 6 

2 8 


54 


34.5/41 


30 


0.94 ±0.01 


31.9 ±2.8 


3.6 ±0.0 




16.9l 3 ;° 


78 


110.5/41 


31 


091 +0.02 


35.9+ 3 ' 2 


3.8 ±0.1 


i3-4±i:i 




43 


45.7/41 


32 


1 06+ 03 


10.6+1;° 


3.7 ±0.0 


12.9 ±0.9 


9 0+ 3 ' 2 
y.u_ 2 g 


48 


71.7/41 



See Tab. ti for an explanation of the fit parameters. The interstellar equivalent column was fixed at Nn = 7.2 X 10 cm 



is a large contribution from the power law in the low-energy 
PHA channels. We therefore expect a strong systematic cor- 
relation between the power law flux, the disk temperature, 
and the disk normalization. In Fig. ^| we display the 99% 
error contours computed for all spectra from the monitoring 
campaign. A strong correlation between the parameters that 
is consistent with the general trend seen in Fig. (?]d is appar- 
ent. This correlation within the individual datasets strongly 
points towards our second argument, that at least part of the 
correlation is due to systematic errors in our model. We note, 
however, that the total variation of the parameters is much 
larger than the uncertainty of the individual data points. It 
is possible, therefore, that at least part of the fcTj n -Adi s k cor- 
relation is indeed real. Observations with instruments with 
a lower energy boundary, however, are needed to verify this 
correlation. 

Apart from a 50 d periodicity which is seen in many 
weak ASM sources (Benlloch-Garcfa et al. 2000, in prepara- 



tion) and is likely due to the nodal precession of the RXTE 
orbit, no significant periodicity was found in the ASM data. 
The overall O(10%) flux variability exhibited on the long 
time scales of the ASM and our monitoring observations 
is comparable to the variability observed on time scales of 
^ 1 ksec, as discussed in paper I. In paper I we attributed 
the variability of LMC X-l to accretion rate variations in a 
wind-fed accretion system. 



5 THE SUPERORBITAL PERIOD OF LMC X-3 

The 100-250 day periodicities seen in the ASM lightcurve 
of LMC X-3 are significantly longer than the measured or- 
bital period of 1.7 days. For the case of Her X-l, the 35 day 
super-orbital period (Giacconi et al., 1973 ) was interpreted 
almost immediately as the result of periodic obscuration by 
a precessing accretion disk which is tilted with respect to the 
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Figure 7. Correlation between the spectral parameters for LMC X-l. Note that the scales of the plots differ from those of Fig. 



plane of the binary system (Katz 1973). The large amplitude 
of the X-ray flux (a factor of 30 between the ASM count 
rate during the "main on" and the "low state" is typical; 
see, e.g., Coburn et al., 1999), as well as numerous other ob- 
servations (see Scott & Leahy 1999 and references therein), 
imply a near edge-on (i ~ 85°) disk. In addition, a frac- 
tion of the central X-ray flux of Her X-l is likely scattered 
into our line of sight by an extended corona (Stelzer et al. 
1999 and references therein). As discussed in the introduc- 
tion, possible driving mechanisms for the disk-warp include 
r adiation pres s ure due to the luminous central X-ray source 
( IPringle 1996|; Maloney et al. 1996| ; |Maloney fc Begelman 



1997 



l^alonejy et al 



1998 



torques exerted by an accretion 
disk wipd (ftchandl kr Meyer 1 994;_ ^chandl 199ri), or tidal 



forces ( Larwood 1998. and references therein). 



Fo r the case of LMC X-3, if the super-orbital variabil- 



ity wer ; also due to changes in the accretion disk geome- 
try, we should be able to observe strong variations in JVh 
as are observed for the turn-on of the Main- On state of 



Her X-l (Kuster et al. 1999). Were the large intensity drops 



in LMC X-3 by a factor of more than five at 5 keV seen 
solely due to variations in Njj_, we would expect columns on 
the order of Nh > 5 x 10 23 cm -2 blocking our line of sight. 
Although such columns are readily detectable with RXTE, 
they were not observed. Given the generally good model fits 
to the observed data and satisfactory residuals, only small 
changes of Nh above the interstellar value appear to be al- 
lowed by the data. We conclude that the strong intensity 
variations of LMC X-3 are not due to the accretion disk 
blocking our line of sight. 

LMC X-3 is spectrally similar to the X-ray binary 
V1408 Aql (=4U1957+11), which also shows evidence of 
long term (100-300 day) periodicities (Nowak & Wilms 



199S). Coupled with evid ence from recent optical observa- 
tions ( Hakala et al. 19~99] ), we have argued that V1408 Aql 



also exhibits evidence of a warped processing disk ( Nowak & 



Wilms 1999). The ASM flux of V1408 Aql, however, varies 



only by a factor of ^ 2 (consistent with warp inclination 
angle changes of ±20°; Nowak & Wilms 1999), which is far 
less than the factor of <; 5 ASM flux variations seen for 
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Figure 8. Same plot as Fig. |?jd, but including 99% error contours 
for each of the fits. The error contours indicate a strong systematic 
correlation between kT ln and Ajisk- 



LMC X-3. Such large amplitude variations would be very 
difficult to achieve unless the inclination of the LMC X-3 
system were comparable to that of Her X-l'f. Given the lack 
of any strongly detected variability on the orbital time scale 
in the ASM lightcurvo of LMC X 3> thie eoome unlikolyi 



Asld 



ism w«pr 



spprtrngmpy yiplHg strong inrtirp- 



tions that LMC X-3 is periodically transitting into and out 
of the canonical low/hard X-ray spectral state of galactic 
black holes. Such transitions are typical for black holes as 
they drop below ~ 5% of their Eddington luminosity (Nowak 
1995; and references therein). If the peak ASM flux corre- 
sponds to ~ 30% LEdd, then the observed transition flux 
threshhold of ~ 0.6 cps in the ASM is roughly consistent 
with 5% of the Eddington luminosity. The question then 
arises of why LMC X-3 apparently shows evidence for large 
accretion rate variations, whereas LMC X-l does not and 
V1408 Aql might not§. 

One possible explanation for the periodic variability 
seen in LMC X-3 is a Compton heated accretion disk w ind- 
driven limit cycle, as discussed by Shields et al. ( 1986} ). At 
large radii, in the optically thin upper atmosphere of the 
outer disk, the central X-ray source will heat the disk tem- 
perature to the inverse Compton temperature given by 



AkTic = I F u hv dvj / F v dv 



(2) 



where F v is the central source flux per unit frequency ( Begel- 



Van cer Klis et al. (h983), however, argue for an inclination of 



t 

LMC X-3 comparable to that which we inferred for V1408 Aql, 
i.e., i ~ 70°. Recent X-ray observations of LMC X-3 (Boyd 1999, 
priv. comm.) show tentative evidence of an orbital modulation. 
This might indicate some warping of at least the outer disk edge 
(see, for example, Dubus et al. 1999), in conjunction with the 
other long term variability mechanisms discussed below. 
§ V1408 Aql historically has had far fewer X-ray observations 
than LMC X-3; therefore, the range of its X-ray spectral varia- 
tions is far less certain. 



man et al. 1983; |Begelman fc McKee 1983| ; [Shields et al 
1986) ; [Schandl fc Meyer 1994| ). If the sound speed for the in- 



verse Compton temperature exceeds the escape velocity in 
the outer disk, a coronal wind can be driven off of the disk. 
For this to occur, the outer disk must have a radius 



Rout ^ 10 1 cm 



M 

V fcTic ) { M 



/3keV\ 



(3) 



LMC X-l, which is likely a wind-fed accretion system (pa- 
per I, and references therein), probably has a much smaller 
circularization radius than this value, perhaps substantially 
so, and therefore is not expected to be driving off a wind 
from its outer disk. V1408 Aql, assuming a 2 Mq primary 
and a 1 Mq companion, has a disk circularization radius of 
~4x 10 10 cm, and also is not expected to be driving off a 
substantial wind (although it may show signs of a photoion- 
ized atmosphere; Nowak & Wilms 1999). For the measured 
system parameters of LMC X-3, the disk circularization ra- 
dius is ~ 1.6 x 10 11 cm, and therefore may be driving a 
substantial Compton-heated wind. 

This opens up the possibility that the central X-ray 
source of LMC X-3 drives off a sufficiently strong wind such 
that it "starves" itself of fuel. The mass accretion rate deficit 
then will propagate inwards on a viscous time scale. The 
central source X-ray flux will decrease thereby shutting off 
the wind, and a limit-cycle can then occur on a viscous 
time scale. This time scale can be of order of 100 days if 
q, the usu al disk viscosity parameter, is ~ 10~ 2 (Shield, 



et al. 198€). One requirement of the instability discussed by 



Shields et al. is that at the peak mass flux of the wind, a 
large fraction, <; 90%, of the total accretion flow is expelled 
from the system. This is consistent with the large amplitude 
variability exhibited by LMC X-3 in the ASM lightcurve. As 
one expects the optical flux to come from large radii and 
the accretion rate fluctuations to propagate inwards from 
the outer disk, this picture is also consis tent w ith the 20 day 
optical lead observed by Cowley et al. (1991). 

T o verify a picture such as discussed by Shields et al. 
([l98§), one must verify that the low/hard flux periods of 



LMC X-3 truly correspond to a 'canonical' low/hard state. 
Additionally, one needs to search for signs of the X-ray 
heated wind. As regards the former possibility, sufficiently 
long RXTE observations during the hard state (> lOOksec) 
should be able to detect the characteristically large ampli- 
tude variability of this state. Furthermore, it recently has 
been noted that low/hard states of galactic black holes, 
as opposed to high/soft states, typically show strong ra- 
dio emission (Fender et al. 1999, Corbel et al. 1998, Zhang 
et al. 1997 and references therein). A previous radio sur- 
vey of LMC X-l and LMC X-3 at 3.5 and 6.3 cm did not 
result in the detection of these objects. The upper limits 
for the flux is 1.5 mjy for LMC X-l. For LMC X-3, an up- 
per limit of 0. 12 mjy at 3.5 cm a nd 0.18 mjy at 6.3 cm has 
been reported (Fender et al. 1998). These upper limits, how- 
ever, were measured during the soft state of LMC X-3. The 
above flux limits do not apply during the recurring episodes 
of the hard state. Scaling from Cyg X-l, LMC X-3 should 
be observable at a flux of ~ 0.05 mjy with a sensitive radio 
telescope during these phases. 
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6 SUMMARY 



ACKNOWLEDGEMENTS 



In this work, we have presented the analysis of an ap- 
proximately two year long RXTE monitoring campaign of 
LMC X-3 and LMC X-l. The primary results of these anal- 
yses are as follows. 

• Both LMC X-l and LMC X-3 can usually be well- 
described by a phenomenological model consisting of a disk 
blackbody with kT n ~ 1 keV and a soft (F > 3) power law 
extending down to the lowest PHA channels. 

• LMC X-l shows 2-20 keV energy flux variations with 
root mean square (rms) variability of 0(10%) on time scales 
of weeks to months. This is comparable to the rms variability 
exhibited on time scales ^ 1 ksec, as discussed in paper I. 
As elaborated upon in paper I, we attributed the variability 
to accretion rate variations in a wind-fed accretion system. 

• Accompanying the flux variations in LMC X-l are cor- 
relations among the phenomenological disk blackbody tem- 
perature, disk blackbody normalization, and power law flux. 
These correlations could be due to variations near a soft 
state/hard state transition region. A strong possibility, how- 
ever, is that these correlations are due to a systematic de- 
pendence of the phenomenological fit components upon one 
another. 

• LMC X-3 shows far larger 2-20 keV flux variations 
than LMC X-l, i.e., a factor of 0(10) luminosity variations 
on 100-200 day time scales. These fluctuations are quasi- 
periodic, albeit with time-varying periods. At the begin- 
ning of the timespan covered by the ASM lightcurve, ~ 100 
and ~ 180 day periods are present, whereas at the end of 
the timespan covered by the ASM lightcurve, ~ 160 and 
~ 240 day periods are present. 

• The flux variations in LMC X-3 are correlated with 
clear spectral variations. At high ASM count rates, the 
flux variations are predominantly associated with varia- 
tions of the phenomenological disk blackbody temperature, 
kTi n , while maintaining a relatively constant disk blackbody 
normalization, Adiak- As the ASM count rate drops below 
« 0.6 cps, the disk blackbody temperature drops, the disk 
blackbody normalization increases, and the power law pho- 
ton index hardens to T ~ 1.8. 

• The spectral variations in LMC X-3 are consistent with 
recurring state transitions between the canonical high/soft 
and low/hard X-ray states of galactic black holes. One possi- 
ble mechanism for explaining the observed periodicities and 
spectral variation s is a wind driven limit-cycle, as discussed 
by Shields et al. fll986| ). 



As discussed above, there are several ways of determin- 
ing whether the low flux/hard spectral states do represent 
canonical low/hard states of galactic black holes. First, a 
relatively long (> 100 ksec) RXTE observation can deter- 
mine whether the fast time scale variability is that expected. 
Second, one can search for an increase of the radio flux con- 
current with the transition to the low/hard state, as has 
been seen for GX 339-4 (Fender et al., 1999) and Cyg X-l 
(Zhang et al., 1997). Evidence of a Compton-heated wind 
should also be searched for in LMC X-3. Sufficiently long 
observations with a sensitive X-ray spectroscopy instrument 
may reveal spectral evidence for a wind (e.g., emission lines 
in the 1 keV region) that may be driving a limit cycle in the 
LMC X-3 system. 
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